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We describe the structure and dynamics of turbulence by the scale-dependent perceived velocity 
gradient tensor as supported by following four tracers, i.e. fluid particles, that initially form a 
regular tetrahedron. We report results from experiments in a von Karman swirling water flow and 
from numerical simulations of the incompressible Navier-Stokes equation. We analyze the statistics 
and the dynamics of the perceived rate of strain tensor and vorticity for initially regular tetrahedron 
of size tq from the dissipative to the integral scale. Just as for the true velocity gradient, at any 
instant, the perceived vorticity is also preferentially aligned with the intermediate eigenvector of 
the perceived rate of strain. However, in the perceived rate of strain eigenframe fixed at a given 
time t = 0, the perceived vorticity evolves in time such as to align with the strongest eigendirection 
at t = 0. This also applies to the true velocity gradient. The experimental data at the higher 
Reynolds number suggests the existence of a self-similar regime in the inertial range. In particular, 
the dynamics of alignment of the perceived vorticity and strain can be rescalcd by to, the turbulence 
time scale of the flow when the scale rg is in the inertial range. For smaller Reynolds numbers we 
found the dynamics to be scale dependent. 



I. INTRODUCTION 

The fundamental question of how small scales are generated in turbulent flows is not fully understood^. 
From a mathematical point of view, small scales are produced by the temporal amplification of the velocity 
gradients when following a fluid clcmcnt»2r™. While the velocity gradient tensor is in principle available in 
direct numerical simulations (DNS), the rare, very intense fluctuations of the gradients make it numerically 
very challenging to accurately determine their statistical properties^. In laboratory flows, of course the 
physics guarantees that all scales are resolved. However, measurements of the full velocity gradient tensor 
are in practice extremely difficult^—. 

To address the question of the generation of small scales, one needs to know the evolution of the flow when 
following a fluid clement. The recent development of Lagrangian measurement technology now allows one to 
experimentally investigate this issued. Nevertheless, the precise measurement of the velocity gradient tensor 
at high Reynolds numbers will for the foreseeable future remain immensely challenging due to temporal and 
spatial resolution constraints. 

Because of the difficulty in determining reliably the velocity gradient, a fruitful alternative approach 
consists in investigating a coarse-grained version of the velocity tensor over a region of size r$. In addition 
to the possibility of determining the flow evolution from a Lagrangian point of view, and thus of addressing 
fundamental questions about the generation of small scales, this approach enables the study of the scale 
dependence of flow properties, especially in the inertial range of scales where the dynamics is believed to 
be determined by inertia only, and dissipation can be considered as negligible. Several approaches can be 
used to define a coarse-grained velocity gradient tensor. One method consists in measuring the true coarse- 
grained tensor by averaging the velocity gradients over a well resolved rcgioni^. The alternative strategy 
followed in this work rests on a reduced descriptio n 13 ^ 4 , based on defining a "perceived velocity gradient 
tensor" (denoted from now on by M(£) or Mij{t)) that is supported only by four fluid elements, initially 
separated by a distance r$ from each other, thus forming a regular tetrahedron. The tensor M(i) reduces 
to the true velocity gradient tensor m(t) when tq — > 0, or more precisely, when the scale r$ is much smaller 
than the Kolmogorov scale, r\, the smallest length scale in the flow. One physically important feature of 
this description is that it provides insight not only on M, but also on the geometry, i.e., the shape of the 
tetrahedron, and on its time evolution. As we see later, the two aspects are strongly couple d 13 i 14 . Thus, 
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the approach based on the perceived velocity gradient tensor and its evolution significantly differs from the 
one based on the true coarse-grained velocity gradient^. ^From a theoretical point of view, the description 
used here replaces a continuous field by a relatively small number of degrees of freedom. It is one of the 
hopes of the present approach that such a conceptual simplification could help in devising simple tractable 
models, which would help shed new light on turbulence. 

The velocity gradient tensor has been shown to possess the property that instantaneously vorticity aligns 
with the intermediate stretching direction (eigenvector) of the rate of strain tensor^£~— , whose eigenvalue is 
predominantly positive. This was surprising as vorticity was expected to align with the strongest stretching 
direction. Already Taylor conjectured 1 "^ that the turbulent cascade mechanism relies on the amplification 
of vortices of scale ro by stretching that subsequently leads to breakup of the vortices into smaller ones. Our 
recent results^ 4 - provide new insight on this conjecture. By analyzing the time evolution of tetrahedra of size 
ro, we found that the perceived vorticity tends to align with the earlier direction of the strongest stretching, 
corresponding to the largest eigenvalue of the perceived rate of strain tensor. For ro in the inertial range, 
this evolution depended only on t/t , where to = (^o/ 6 ) 1 ^ 3 i s the time scale of eddies of size r in a turbulent 
flow with energy dissipation rate s. From here on, we use "vorticity" and "rate of strain" both for the 
perceived as well as real quantities, the difference can be deduced from the scales given. 

Here, we analyze the structure and the dynamics of M, and its dependence on scale, thus extending 
our previous investigations^ 4 -. Our experimental results at R\ = 350 demonstrate that the dynamics of the 
alignment of vorticity with the eigenvectors of the rate of strain, measured at the time where the tetrahedron 
is regular, depends on r through the rescaled time t/t n . This suggests the existence of a self-similar regime 
in the inertial range. At smaller Reynods numbers, we found the dynamics to be scale dependent. The 
numerical results show that when ro — > the properties of M extrapolate smoothly to those of the true 
velocity gradient tensor, and the dynamics of alignment between vorticity and strain is qualitatively similar 
with one notable difference: the characteristic time scale for tetrahedra with size ro < r\ is of the order of 
tki the Kolmorogov time scale. 

We begin, Section UH by describing the construction of the velocity gradient tensor M. Our experimental 
and numerical techniques are explained in Section IIHI Section IIVI then presents our results concerning the 
evolution of the shape of the tetrahedra. In Section [Vj we discuss the instantaneous correlations between 
strain and vorticity at any instant, and their dependence on the scale ro. The dynamics of M, in particular 
the alignment between vorticity and the initial eigenvectors of strain are presented in Section I VII which 
includes a subsection of simplified theoretical analysis of the alignment process at short times. Finally, 
Section [VIII contains our concluding remarks. 



II. DEFINITION OF THE PERCEIVED VELOCITY GRADIENT TENSOR M 

Given fluid velocities at four points in a flow, we compute the perceived velocity gradient M as follow o 13 i 14 . 
Denoting the positions of the four points as x a (i), and their velocities as u a (t), a = 1, ...,4, we express the 
positions and velocities with respect to the center of mass as x' a = x a — x°, and u' a = u a — u°, where 
x° = ^(x 1 + x 2 + x 3 + x 4 ) and u° = | (u 1 + u 2 + u 3 + u 4 ). In full generality, My is defined by minimizing 



In the equation defining the quantity to be minimized, K, the indices a (= 1, 4) refer to the index of the 
point, and i, j (= 1, 3) to the component of the vector. By differentiating with respect to My, one finds 
immediately that the condition to minimize K is: 



the quantity: 
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where the matrices g and W arc defined by: 
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Under the condition that g is invertible, the matrix M defined by: 

M = g _1 W, (4) 

provides the best fit approximation of M based on four points. 

As the flow is incompressible, we impose the constraint that tr(M) = 0. This can be done by determining 
M, as explained in the previous paragraph, and subtracting gtr(M)5y-. Alternatively, returning to the 
minimization condition, Eq. ([!}, the trace can be imposed by adding a Lagrange multiplier, £,^2 i Mu to 
Eq. ([1]). The two definitions, although technically slightly different, lead to equal statistical properties, and 
to identical physical conclusions. 

The construction used here, based on Eq. (j4]), explicitly requires the matrix g to be invertible. As a 
symmetric matrix, g is diagonalisable. Its eigenvalues gi are positive, and yfgl characterize the extent of the 
tetrahedra in the i th eigendirection of g. The matrix g becomes non-invcrtiblc when one of its eigenvalue 
tends to zero. Thus, the condition that g is invertible means physically that the four points of the tetrahedron 
are not coplanar. The technical difficulty associated with the minimization of the quantity K for highly 
flattened tetrahedra can be solved using standard Singular Value Decomposition algorithms 2 ^. 

In this work, we followed tetrahedra that are initially regular (close to regular in experiments), so as to 
identify flow properties at a single scale. The alignment dynamics studied here occurs before the tetrahedra 
become strongly deforme d 21 i 22 or we did not use the highly deformed tetrahedra in the analysis. As shown 
in Section ITVl the fraction of highly-flattened tetrahedra is negligible up to io/4 and only approximately 10% 
are highly-deformed at £o/2- As will be shown in later chapters, the most interesting dynamics occur before 
to/2. To further reduce the effect of tetrahedron deformation on the accuracy of determination of M, we 
excluded highly deformed tetrahedra from our statistics of M (see details in Section IIV|) . This comes at the 
cost that the number of samples in our statistics is slightly reduced. 

Similar to the treatment of the true velocity gradient, it is convenient to decompose M as the sum of a 
symmetric, strain-like part, S, and an antisymmetric, vorticity-likc part, J~2: 

S = i(M + M T ) ; n = i(M-M T ). (5) 

The symmetric matrix S describes the local straining motion. It is characterized by three real eigenvalues, 
Ai, A2 and A3, associated with three eigenvectors e%. In the following, the three eigenvalues are sorted in 
decreasing order: Ai > A2 > A3. The antisymmetric matrix, f2, describes the local rotation: Qy = iey-feWfe, 
where cijk is the completely antisymmetric tensor, so that H • x = ^u> x x. In the following, we characterize 
to by its norm \u>\ and its direction e w , with |e w | = 1. In the limit of very small tetrahedra (tq <C 77), the 
tensor M reduces to the true velocity gradient tensor m, defined by my = djU iy and the vector ui defined 
above reduces to the usual definition of vorticity, oj = V x u. 

III. EXPERIMENTS AND NUMERICAL SIMULATIONS 
A. Experiments 

Using image-based optical particle tracking, we measured experimentally the trajectories of tracer particles 
seeded in a swirling water flow between two counter-rotating baffled disk o 23 ' 24 . Measurements were done 
in the center of the apparatus, where the residence time of particles is the largest. The maximum Taylor- 

microscale Reynolds number is R\ tn 10 3 , where R\ = \J\hu'^ jve with u' = W '(it 2 + it 2 + u 2 ) /3 being the 

fluctuation velocity and v the kinematic viscosity. The symbol () refers to the average of the fluctuating 
quantities in the flow. Measurements were carried out by using three high-speed CMOS cameras and high 
repetition rate frequency doubled Nd:YAG lasers. Our particle tracking algorithm 2 ^ allowed us to follow 
simultaneously hundreds of particles in this intense turbulent flow. The measured trajectories, however, 
were frequently interrupted due to experimental artifacts such as the fluctuation of the illumination light 
intensity, the background and electronic noise. For multi-particle measurements, such as tetrahedra, it is 
very important to connect these interrupted trajectory segments in order to obtain better statistics at long- 
times. We therefore re-connected the trajectories in position-velocity spaced, which were then smoothed 
and differentiated to obtain instantaneous particle velocities^. 

In this paper, we focus on results from a measurement at R\ — 350, at which the scale separation between 
the integral scale L = u' 3 /e and the dissipative scale r\ = (i> 3 /e) 1 ^ 4 is L/rj ~ 800. As shown previously 2 ^, 
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the temporal and spatial resolutions are sufficient to accurately measure the second moments of acceleration 
statistics. 

^From experimentally measured particle trajectories, we first identified instantaneous configurations of 
four points, separated from each other by a nominal distance fo, within a tolerance of ±10%. Once identified, 
the trajectories of the tetrahedra were followed for as long as all four points remained in the field of view. 
Due to the limitation of particle seeding density in experiments, the data reported are mainly for tetrahedra 
with initial sizes between L/17 — L/5, or 50 — I8O77. In this range of scales, we followed 10 7 — 10 s tetrahedra. 
The small tolerance of inter-particle distances at the nominal values results in initially nearly isotropic 
tetrahedra— . 

The particle trajectories used in this work were measured in an observation volume of size approximately 
(2.5 cm) 3 located at the center of the apparatus, where the flow has been documented to be very close to 
isotropic*— Deviations from isotropy may be estimated from the well-known second order velocity structure 
functions, defined as D LL (r) = ([(u(x + r) - u(x)) • (r/r)] 2 ) and D NN (r) = {([(u(x + r) -u(x))] 2 ) -D LL }/2, 
where r = |r| is the separation distance. For isotropic turbulence, £>ll and I? at at depend only on r and they 
satisfy the following relations (see e.g.,—): Dnn{t) / Dll{t) = 2 for r <C rj and Dnn{t)I 'Dll{t) = 4/3. For 
the von Karman flow in our experiment, due to the axisymmctric geometry, the statistics could depend on 
the angle 8 between the axis of rotation and the separation vector r. We found that for 40ry < r < L/3, the 

ratio (3/4)Djvjv( , r, 6)/Dz,L(r, 6) varies as a function of 9 in the range 0.8 1.3, but averages out to be very 

close to 1, the isotropic value, when considering all possible orientations of the separation. This provides 
evidence both that the dependence on the angle 6 is weak, and that isotropy is recovered when averaging 
out over all possible direction s 30 ' 31 . Further evidence that the anisotropics in our flow are not a concern 
comes from the fact that most experimental statistics derived from tetrahedra are in good agreement with 
DNS results using perfectly isotropic tetrahedra, see also— The small deviation from precise isotropy in 
initial shapes, however, did give rise to discrepancies for certain quantities sensitive to initial shape, which 
we will discuss in detail later. 



B. Direct Numerical Simulations 



We used a spectral code to solve the Navier-Stokes equations for the (Eulerian) velocity field u(x, i): 

9 t u(x, t) + (u ■ V)u(x, t) = - Vp(x, t) + ^V 2 u(x, t) (6) 
V-u(x,t)=0 (7) 

Eqs. ([6]) and ||7J) are integrated in a periodic box (size 27r), with up to 384 3 modes. Energy is injected into 
the flow at large scale by letting the low wavenumber modes |fc| < K (here, K = 1.5) evolve according to the 
Euler equations truncated to the shell \k\ < K . This maintains a constant amount of energy at the largest 
scales— Because the simulated flow is highly isotropic, we define the turbulence parameters based on the 
x component of velocity. The Reynolds number is 

Rx - ^ (8) 



where the Taylor micro-scale defined as 



9 ((d x u x )2)- 



(9) 



(it 2 ) 1 / 2 . The Kolmogorov length scale is rj = (z^/e) 1 / 4 , where e = v(djUidjUi) . The 
available resolution has allowed us to simulate reliably turbulent flows up to a Reynolds number R\ = 170. 
Quantitatively, the largest wavenumber faithfully simulated, k max , is such that the product k max x 77 w 1.4. 

Defining the integral length scale by L = u' 3 /e, the ratio L/rj is approximately 300 at R\ = 170. We note 
that L defined here is roughly twice as large as the correlation length scale of the Eulerian velocity field: 



Lcorr - 



— / k~ l E{k)dk w L/2, (10) 
u Jo 



where E(k) is the energy spectrum. 
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With the Eulerian velocity field determined by solving numerically Eqs. ^ and (J7J), we also followed the 
motion of tetrahedra, whose vertices are tracer fluid particles, evolving according to: 

^ = u(x(t),t). (11) 

To this end, the values of u, known after each time step on a regular grid of collocation points, are interpolated 
at the location x(t) of the Lagrangian tracers, using accurate third order schemes^ 3 .. 

To study the dynamics using tetrahedra, we initialized tetrahedra of particles mutually separated by a 
distance ro, which we varied over the full range of available scales, r\ < ro < L. The statistical properties at 
a given scale r have been determined by following at the minimum 125, 000 tetrahedra. 



IV. DEFORMATION OF THE TETRAHEDRA 



We consider here the deformation occuring over a short time scale, while the process of alignment between 
the perceived vorticity and strain is taking place. To quantify the deformation of tetrahedra, it is convenient 
to introduce the tensor g, defined in Eq. ((3]), and its eigenvalue o 14 ' 21 ! 22 ' 34 ! 35 . The trace of g, which is the 
square of the radius of gyration of the four points, quantifies the overall size of the object. The tensor 
I = g/tr(g), whose trace is 1, characterizes the shape of the object. Specifically, because I is a symmetric 
matrix, it can be diagonalized in an orthogonal basis formed by its eigenvectors en. In the following, we 
arrange the eigenvalues of I in decreasing order: Ii > I<i > I3. Since I\ + L 2 + I3 = 1, a regular tetrahedron 
corresponds to I\ = I2 = I3 = 1/3, while a nearly co-planar configuration gives I3 « 0. 

The matrix g is related to the moment of inertia tensor, J, familiar in the context of mechanics of solid 
bodies, byi^: 

Jij =tT(g)Sij -gij. (12) 



A. Alignment of tetrahedron geometry with the eigenvectors of the rate of strain 



The dynamics of tetrahedra is strongly influenced by the local rate of strain. In particular, it has been 
demonstrated^ that the stretching in the direction of the largest eigenvalue of the strain, ei, leads to an 
elongation of the tetrahedron in this direction, and, by conservation of angular momentum, to an alignment 
of the vorticity u(t) in the direction ei(0). The results presented here generalize this observation: at short 
times, the deformation of the tetrahedra is very well aligned with the eigendirection of the eigenvectors of 
the rate of strain tensor. 

Fig. [T] shows how the principal axes of initially regular tetrahedra, en(t), align with the eigenvectors of 
the rate of strain, e;(0). For initially isotropic tetrahedra, as in DNS, the matrix I is diagonal at t = 0. A 
very small deformation of the tetrahedra due to the strain breaks the isotropy, and leads to the immediate 
alignment of the principal axes of the tetrahedra with the strain eigenvectors. This is reflected in Fig. [1] 
by the perfect alignment between efc(0) and eik(t) at very small times t (t = + ). By t = + , we refer 
to a strictly positive, but very small value of time. The tetrahedra studied experimentally are not exactly 
isotropic initially. Instead, the eigenvalues of I are slightly different from each other: 7j = i + 5Ii, where 
\6Ii\ < 1. At t = 0, the corresponding three eigenvectors of I are uncorrclatcd with the strain eigenvectors, 
i.e., ([ejfc(0) • e 3 (0)] 2 ) « 1/3. This explains the difference at short times seen in Fig. [1] between the curves 
obtained from DNS and from experiment. On the other hand, under the action of the strain, the tetrahedra 
deform quickly and their principal axes align with the eigenvectors of the strain in a time of approximately 
to/10, after which the tetrahedron evolution is independent of initial shape. We note that over the range 
of scales shown here, L/16 < ro < L/2, the dynamics of alignment is essentially self-similar, once expressed 
in units of to- In a time of order to, the axes eik{t) lose any memory of the initial alignment with ej(0); 
all values of ([ejk(t) ■ ej(0)] 2 ) relax towards 1/3. This is consistent with the fact that to, which is the time 
introduced in Kolmogorov theory, is the proper correlation times of an object initially with a scale ro. Small 
deviations are visible for DNS data of ro = L/16 sa 20i], which is known to be the lower limit of the inertial 
range. 
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FIG. 1. (Color online) Deformation of tetrahedra: alignment of the principal axes of the tetrahedra with the 
eigenvectors of the strain. The ensemble averages {[en(t) ■ ej(0)] 2 ) as a function of time for (a) i = 1 (longest axis), 
(b) i = 2 (intermediate axis) and (c) i = 3 (shortest axis of the tetrahedron). As expected, the i th axis of the 
tetrahedron, en(t), aligns perfectly at very short times with the i th eigenvector of the strain, ei(0). After a time of 
order ~ to, the alignment relaxes and the axes of the tetrahedron do not show any particular alignment with any of 
the eigenvectors of the strain. The DNS data correspond to R\ — 170, whereas the experiments to Rx = 350. 



B. Flattening of the tetrahedra 

Along with the stretching of the tetrahedra in the direction ei(0), documented in particular ir>ii, one 
expects a strong compression in the direction 63(0), corresponding to the most negative eigenvalue of the 
perceived strain. This has been shown to lead to very significantly flattened tetrad s 21 ' 22 ' 35 . In particular, the 
probability distribution function (PDF) of I3 is sharply peaked around at late times, t/to > 1. As explained 
in Section [XXI this can be the cause of a serious limitation in our construction of the perceived velocity 
gradient, as implied by Eq. A good understanding of the dynamics leading to flattened tetrahedra is 
thus of interest not only for its own sake, but also in relation to the dynamics of alignment between perceived 
vorticity and strain. This subsection therefore discusses the flattening of tetrahedra, and demonstrates that 
the formation of very flat configurations is not faster than the alignment process between vorticity and 
strain. 

To this end, Fig. [2] shows the evolution of the PDF of I3 (a), and the evolution of the probability that the 
value of I3 is less than a very small number, ^3 < 10~ 3 (b). Initially, the probability of ^3 is peaked near 
^3 = 1/3. In DNS, which starts with precisely isotropic tetrahedra, the PDF is originally a 5-function at 
J3 = 1/3. In experiments, where we select particles that form nearly isotropic tetrahedra, the initial PDF 
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FIG. 2. (Color online) Flattening of tetrahedra. (a) PDF of 7 3 at several times, t/t = 0.07, 0.13 and 0.21, for 
different initial tetrahedron sizes (ro from L/16 to L/2) in both DNS (R\ = 170) and experiment (R\ = 350). The 
observed peak of the PDF decreases when time increases. At t = 0.21to, a non-zero PDF at 73 ~ is observed. The 
probability that 73 < 10 -3 is shown in (b) for several values of ro, as a function of the normalized time t/to- A steep 
growth of the probability P(73 < 10~ 3 ) is seen for t > to/4. The growth of the probability does not depend much on 
the scale ro over the inertial range: from ro = 7/16 ~ 20r/ in DNS and ro = 7/17 ~ 50r/ in experiments to ro = L/2. 



is not a (5-function, but is still narrowly distributed around 73 sa 0.3, as shown before^i. As time increases, 
the peak of the PDF shifts towards smaller values of 73 (Fig. [5^), reflecting a tendency of the tetrahedra 
to flatten. The evolution of the PDF does not depend very much on the initial size of the tetrahedron, as 
seen when comparing DNS and experimental data ranging from ro = L/16 as 20r/ in DNS to ro = L/2. A 
noticeable fraction of the tetrahedra becomes flattened after i > to/4 (see Fig. Hfb)), consistent with the 
sharply peaked distribution of 73 observed at later time o 21 ' 22 i 35 . ^From a flat configuration, with all points 
in the same plane, it is not possible to extract any information on the variation of the velocity in the plane 
transverse to the plane. As explained in Section |TT] however, this does not prevent us from computing M 
using Singular Value Decomposition. However, it points to a limitation on the accuracy of the procedure. 
Therefore, in the results reported here, we excluded tetrahedra with 73 < 10 -3 from the statistics of M. As 
seen from Fig.[5fb), this leads to at maximum a 10% reduction in the number of samples at time t/to = 0.5, 
while most of the interesting dynamics occurs before this time. 

The results of this section thus show a very strong alignment between the axes of the tetrahedra and 
those of strain, sec Fig. [TJ and that the dynamics of flattening of the tetrahedra happens over a time of the 
order of to/4, sec Fig. [21 In addition, we find that the dynamics is essentially self-similar, insofar as the 
time-dependent correlations presented here are independent of scale, once time has been expressed in units 
of to- This is consistent with earlier results 2 ^ demonstrating that the evolution of the shape factors, 7j(t), 
is also self-similar. The ensuing picture is thus that in a statistical sense, initially regular tetrahedra of size 
ro evolves with a time scale ~ (r^/e) 1 / 3 , when ro in the inertial range. 



V. PROPERTIES OF M: SCALE DEPENDENCE (INSTANTANEOUS STATISTICS) 

Previous investigation o 6 ' 10 ' 16 ! 17 ' 36 of the structure of the true velocity gradient tensor, m, reveal two main 
properties of the instantaneous statistics of m: 

(i) The direction of the vorticity is preferentially aligned with e%, the direction of the intermediate 
eigenvalue of the rate of strain, s = ^(m + m T ), but does not show any particular alignment with ei, the 
strongest eigen-dircction of the rate of strain. 

(ii) The intermediate eigenvalue of the rate of strain, A2, is mostly positive, which implies that the product 
of the eigenvalues of s is negative and thus ensures vortex stretching (cj • s • u)) > 0, an important property 
of turbulent flows^Z-. 

It is natural to ask what are the instantaneous properties of the perceived velocity gradient M. 
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On general grounds, one expects that the properties of M will depend on the size tq. For r < r/, M 
coincides with the true velocity gradient tensor. As the separation ro between the points of the tetrahedron 
increase, the velocities at the four points become less and less correlated. For ro 2> L, M is not expected to 
reflect any particular property of the turbulent dynamics. On the other hand, one may expect that in the 
inertial range, some properties become scale-independent, consistent with the notion that the flow shows 
some degree of self-similarity 



A. Alignment between vorticity and strain 

The PDFs of the cosines of the instantaneous angles between and ej , are shown in Fig. 03 The PDF 
for i = l (Fig. (Ha) and (b)) hardly varies through scale, and remains nearly uniform: P{\e\ ■ e w |) ~ 1, 
consistent with previous observations with the true velocity gradient tensor— In Fig. [3jc) and (e), the 
alignment corresponding to the true velocity gradient tensor m (i.e., ro = 0) are the ones with the largest 
variations. As scale increases, the alignment between e 2 and becomes less pronounced, as shown by the 
observed decrease of the PDF at |e2 ■ e u \ = 1. Similarly, the tendency of e w to be perpendicular to e3 
diminishes when scale increases. As the scale ro approaches L, the PDF becomes flat. This tendency can 
be clearly seen in Fig. 21 which shows the averaged values ((ej -e^) 2 ). When ro increases from dissipative to 
integral scales, ((e u • e 2 ) 2 ) continuously decreases from the larger magnitute for the true velocity gradient 
tensor (shown by the horizontal dashed lines in Fig. 2]) towards 1/3, which corresponds to two randomly 
selected vectors. In the same spirit, ((e w • e3) 2 ) continuously increases from the lower magnitude of the 
velocity gradient tensor to 1/3 when r increases. In Fig. 2] are shown both the numerical data (R\ = 170, 
open symbols) and the experimental data (R\ = 350, filled symbols). Remarkably, in both Figs. [3] and 
21 the lower Reynolds number data show a continuous variation of the statistics of the cosines e 2 • e u and 
e3 • e^ with scale. The variation is particular strong when ro decreases towards the dissipative range, i.e., 
f° r To ^ L/16 ~ 20rj, In comparison, the higher Reynolds number data suggest that over the accessible 
range of scales, the properties of the flow are essentially scale independent. This suggests that at sufficiently 
large Reynolds numbers, the dynamics of M may show some self-similarity, as far as the relation between 
vorticity and strain is concerned. This remains to be studied further with data at larger Reynolds numbers. 

The results summarized in Figs. [3] and H] thus show that the alignment properties of e w and e 2 , or the 
orthogonality between e w and e3, are observable throughout the entire inertial range. This is to be contrasted 
with the lack of alignment between e w and ei, which is observed over all the scales that we have studied. 



B. The intermediate eigenvalues of the rate of strain 

We now turn to the statistics of Aj, the eigenvalues of the rate of strain. Whereas it is clear, from the 
relations Ai > A 2 > A3 and Ai + A 2 + A3 = that Ai > and A3 < 0, the sign of A 2 can be either positive 
or negative. To quantify the relative value of A 2 compared to Ai and A3, it is convenient to introduce the 
dimensionless variable j3, defined byii: 

y/X{ + A 2 + A§ 

This definition ensures that — 1 < j3 < 1. 

The PDF of j3 is shown in Fig. [5^ for several values of ro. At small Reynolds number, the PDFs show 
a clear evolution as a function of the scale ro. In the dissipative range (ro = 0, i.e., for the true velocity 
gradient tensor), a strong bias towards high values of (3 is observed. As r increases, the PDFs show a 
reduced bias towards high values of (3. For the highest Reynolds number experimental flow, the PDF of j3 
changes little throughout the limited range of measured scales. As shown in Fig. [SJd, the mean ((3) remains 
positive throughout the entire inertial range, dropping from (f3) w 0.285 for ro = 0, to (j3) ~ 0.07 for ro ~ L 
for the low Reynolds number data. On the contrary, the value for the higher Reynolds number experimental 
flow, in the range 1/17 < r /L < 1/5 is close to 0.13. Figure [5] suggests that the dependence of (j3) on scale 
is relatively weak in the inertial range L/16 < ro < L/2, in particular at the highest Reynolds number. This 
is consistent with the trend observed in subsection IV AI which points to a scale-independent structure of M 
in the inertial range at sufficiently high Reynolds numbers. 

The results of this section suggest a smooth change of the structure of M throughout the range of scales. 
Whereas no particular alignment between e w and ei is observed, irrespective of the size of the tetrahedron, 
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FIG. 3. (Color online) Scale dependence of the alignment of vorticity and strain (instantaneous statistics): The 
PDFs of |ei ■ e„| from both DNS (panels a,c,e) and experiments (panels b,d,f) for tetrahedra with different sizes ro 
from dissipative (ro = in DNS) to integral {L/2) scales. Panels (a) and (b) are for i = 1; (c) and (d) for i — 2; 
and (e) and (f) for i — 3. The lack of alignment of vorticity with ei is nearly independent of scale. The PDF of 
|e„ -e2| sharply peaks at 1 at very small scales, in particular for the true velocity gradient tensor (ro = 0). The effect 
weakens as the scale increases. Similarly, the maximum at of the PDF of |e„ ■ es| at very small scales becomes 
milder as scale increases. The DNS data correspond to Rx = 170, whereas the experiments to R\ = 350. 
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FIG. 4. (Color online) The mean values {(e„(0) •ei(O)) 2 }. The open symbols are DNS data (R\ = 170) and the filled 
symbols are experimental results (R\ = 350). 



the tendency of e w to align with e 2 (respectively to be perpendicular to 63) is the strongest for vq — > 0, i.e., 
for the true velocity gradient tensor m, and decreases for M with increasing tq. Similarly, the intermediate 
eigenvalue of the rate of strain is comparatively more positive for the true velocity gradient tensor m than 
for M at values of r$ in the inertial range. As the scale ro increases towards the inertial length L, the PDF 
of P becomes more symmetric about (3 = 0, as expected from a random matrix. 

Thus, the perceived velocity gradient tensor M in the inertial range shares many essential properties 
with the true velocity gradient tensor m. Our observations thus suggest a continuity in the behavior from 
the dissipative to the inertial range. Remarkably, the data obtained with the higher Reynolds number 
experimental flow suggest that the structural properties of M become scale independent in the inertial 
scale. This should be confirmed by studying flows at even higher Reynolds numbers. 



VI. PROPERTIES OF M: DYNAMICS (TEMPORAL CORRELATIONS) 

The main observation in our recent studyi^ was the dynamical alignment of vorticity with the leading 
eigenvector of the rate of strain, which we were able to study conveniently by investigating the correlation 
between ei(0) and e w (£), i.e., with a time delay. Hence the lack of alignment between e w (t) and e±(t) 
(Figs.|3fa) and (b)) is a consequence of the de-correlation of the vector e! itself in a time t, which counteracts 
the alignment of e w (t) with ei(0). 

A consequence of the previous section is that it is not a priori clear that the alignment of e w (t) with e; (0) 
(i = 2 or i = 3) is going to be self-similar, i.e., that the angle between e u (t) and e,(0) (i = 2 or i = 3) is 
going to depend only on t/to, as it had been found for the angle between e u (t) and ei(0)^. This is due 
to the fact that the PDFs of |e w (0) • ej(0)|, for i = 2 and i = 3 change with the scale ro for small scale 
separation L/rj. 

A. Averaged alignment of the direction of vorticity with the eigenvectors of the rate of strain 

Fig. |5{a) shows that, as time increases, the vector e w (i) tends to become better aligned with e^O). As 
reported previously^, the statistics characterizing the alignment of the two vectors is essentially independent 
of the scale ro, once time t is rescaled by to- The characteristic time scale for this process, Fig. Ela) is of 
order to/5. For the DNS at R\ = 170, the alignment of e u (t) with e2(0), and to a lesser extend with 63(0), 
shows a significant dependence on scale, see Fig. EJb) and (c). Note that for t = 0, the difference was shown 
as a function of ro/L in Fig.(¥] In comparison, the experimental data obtained at R\ = 350 shows a much 
reduced variation as a function of scale (see also Fig. @|. We also notice that ([e 2 (0) • e tJ (t)] 2 ) is almost 
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FIG. 5. (Color online) Dependence of the intermediate eigenvalue of strain on scale ro. (a) PDF of j3, as defined in 
Eq. 1131 for several values of ro. (b) The average value of /3 as a function of scale, in which the dashed line indicates 
(/3) for the true velocity gradient m as ro — > 0. The DNS data correspond to R\ — 170, whereas the experiments to 
Rx = 350. 



constant, i.e., the slope at t = is much smaller than the slopes for ei(0) and 63(0). These features can be 
understood by the elementary considerations presented in subsection IVI Bl 

We note that the time scale of the alignment of e w (i) with ei(0) and 63(0) are comparable to the time 
scale characterizing the deformation of the tetrahedra (see Section HV|) . This is a strong indication that the 
deformation is an important part of the dynamical process, see also 14 ' 38 . A thcorical description of the 
dynamics of M has therefore to take into account in an essential way the deformation of the tetrahedrai 3 -. 

Not only the mean values of [ei(0) • e w (i)] 2 evolve as a function of t/to in a self-similar way, but also the 
PDFs of |ei(0) ■ e w (i)| (sec Figures [TJa), (d), and (g)). This situation is to be contrasted with the evolution 
of the alignment between e w (i) and e2(0) or e3(0). FigurcslTJb), (e), and (h) show that the PDF of the angle 
between e 2 (0) and e w (t) is essentially independent of t/t for < t/t < 0.2, consistent with Fig. [B^b). The 
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(b) and i = 3 (c) as a function of the time delay t, normalized by to- The initial tetrahedron size, ro, are all in the 
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FIG. 7. (Color online) Evolution of the angle between ei(0) and e w (t). The PDF of |ej(0) ■ e u (i)| is shown for times 
between < t/to < 0.2, and at three values of ro: ro = L/2 in DNS (R\ = 170) (a)-(c), ro w L/9 and ro ~ L/17 in 
experiment (R\ = 350) (d)-(f), and ro = L/16 in DNS (g)-(i). Panels (a), (d) and (g) are for i = 1. Consistent with 
previous results^, the evolution of these PDFs is essentially self-similar for ro in the inertial range. Panels (b), (e) 
and (h) are for i — 2. The PDFs hardly change over the range of times shown; however, they depend significantly 
on the value of ro. Panels (c), (f) and (i) are for i = 3. During the period of time shown, the PDFs evolve to peak 
at 0, i.e., e u (t) tends to become perpendicular to e3(0). There is also a moderate dependence on scale ro. 



direction e w (t) is strongly aligned with e2(0) at small values of ro, much less so as the value of ro increases 
towards L. Last, as shown in Figs. EJc), (f), and (i), e u (t) becomes more perpendicular to es(0) within 
t/to ~5 0.2 and the tendency is stronger at small values of ro. 

In summary, the main observation of this subsection is that, contrary to the alignment of e u (t) with ei(0), 
which evolves in an essentially self-similar manner, the dynamics of the alignment of e w (i) with e2(0) and 
e3(0) proceeds in a way that depends on the size of the tetrahedra, at least for our Reynolds numbers, a 
result that was expected in view of the observation that the average values of (ej(0) • e w (0)) 2 show a clear 
dependence on scale. 
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B. Elementary modeling considerations 

The analysis presented below, which aims only at capturing some of the main features of the alignment 
process, rests on several important assumptions. We restrict ourselves here to very short time scales, which 
effectively allows us to neglect the coupling between shape and vorticity. We are merely investigating here 
the evolution of vorticity, at the lowest order in terms of a Taylor-series expansion. Effectively, coupling 
terms between shape and vorticity, as well as feedback between vorticity and strain all would be seen only 
at higher order, independently of the form of the model. The limitation of this approach is that it does not 
allow us to properly describe the coupling between shape deformation and the evolution of vorticity. To 
describe the evolultion of M, models based on the Restricted Euler approximation^— and its subsequent 
elaborations, reviewed hJ£, have been used very often. The resulting equation for the evolution of the 
vorticity is: 



dt 



S. 



= 0. 



(14) 



which can also be obtained directly from the Euler equation. A full description of M requires also an equation 
for the evolution of the rate of strain, S, which involves the second derivatives of pressure (pressure hessian), 
a quantity notoriously difficult to model. While it may be possible to extract information concerning the 
strain evolution from the available data, we restrict ourselves here to the evolution of vorticity, and focus 
exclusively on Eq. [TU 

An evolution equation for the direction of vorticity e w can be readily obtained by taking the double vector 
product of Eq. [TJ] with e w , which immediately leads to: 



^ = S • e w - (e u ■ [S ■ e u ])e u 



(15) 



To proceed, we project e w on the three eigenvectors of S: 



(16) 



The constraint |e w | 2 = 1 implies that Yl^=i c f = 1- Eq. [15] can be simply written in terms of the three 
components Cj. As we are primarily interested here in ((e w • e^) 2 ), we use Eq. [15]for cf: 



dt 



= 2c 



(17) 



The initial dynamics of alignment is best investigated by averaging over many realizations, and by sepa- 
rating the three components in Eq. (I17[) : 



dcf_ 
dt 



= 2(c 2 c 2 (A 1 -A 2 ))+2( c 2 c 2 (A 1 -A 3 )) 



t=o 



(18) 



def 
dt 



t=o 



(cf^(A 1 -A 2 ))+2(c 2 c 2 (A 2 -A 3 )) 



(19) 



and 



d4 

dt 



= -2<c?c§(Ai - Aa)> - 2<cfc§(A a - A3)} 



(20) 



dc 

The sign of (-^-)\t=o is determined solely by the eigenvalues of the rate of strain, with Ai > A 2 > A 3 . We 

dc 2 d 2 

see immediately that (-§*)|t=o > and (-%f-)\t=o < 0. Thus initially vorticity tends to align with ei(0) 

d 2 

and to become perpendicular to e 3 (0). By the following argument we also see that (-^f-)\t=o is smaller in 
magnitude than the other two components. The first term on the RHS of Eq. [19] is negative and has the 
same magnitude as the first term on the RHS of Eq. II 81 hence this term is negative and smaller in magnitude 
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than (-^)\t=o- Also, the second term on the RHS of Eq. [12] is positive and has the same magnitude as the 
second term in the RHS of Eq. [2D] therefore this term is smaller than the magnitude of (-3#)|t=o- Thus 

d 2 

(-j#)|t=o is smaller in magnitude than the other two terms and consequently its evolution is the slowest. 
This agrees qualitatively with the observations reported before. 

The analysis presented above is merely based on the action of strain on the vorticity, which is common to 
essentially all models aimed at describing the evolution of M. Comparing more quantitatively the results of 
the present work with the prediction at short times of more sophisticated models, such as the tetrahedron 
modeli 3 - or the Lagrangian stochastic mode l 42 ' 43 , requires the determination of several constants in the 
models, which is beyond the scope of this work. We nevertheless note that in the specific case of the 
tetrahedron modeU 3 -, the equation of evolution of vorticity reads: 

^ - (1 - a)SijU>i = (21) 

where a, which describes the reduction of nonlinearity, is a parameter that a-priori depends on the scale 
7'o • Our observation that the evolution of M may be self-similar at large Reynolds numbers throughout the 
inertial range suggests that the parameter of the model are independent of the scale ro . 

As explained already, an extension of this calculation to later times requires both the knowledge of the 
interplay between vorticity and the rate of strai n 36 ' 44 ' 45 , as well as the geometry of the tetrahedra. The very 
good superposition of the correlations ((e u (f) • ei(0)) 2 ) in Fig. [6] corresponding to tetrahedra of different 
sizes, while ro is in the inertial range, is consistent with a description of the dynamics of M by such a model 
based on the Restricted Euler approximation or one of its gcncralization o 18 ' 39 ' 41 ' 42 . Indeed, on general 
grounds, M is expected to behave with scale as M oc (e/r 2 ,) 1 / 3 , so the eigenvalues of the rate of strain 
\ scale as A^ oc (e/r 2 ,) 1 / 3 , which in turn is consistent that the dynamics of alignment proceeds with the 
characteristic time scale to, as observed before^ 4 - and in this work. More generally, a dynamical equation of 
the form cx M 2 immediately suggests that the characteristic time scale of the dynamics should be to, 
thus suggesting that all the dynamical process should depend on t/to, at least when ro is in the inertial range. 
While our results on the alignment between e u (t) and ei(0) strongly suggests some universal behavior as 
a function of scale ro , it is not settled by the present work whether the dynamics of alignment of vorticity 
with the other two eigendircctions of strain is also self-similar over a range of scales in the inertial scales. 
It would be very interesting to determine whether a dependence on scale ro of the dynamics of M persists, 
even at very large Reynolds numbers. These issues are important in many aspects for the modeling of the 
flow at scales r , in particular concerning the applicability of models based on Restricted Euler, and its 
generalizations. 



C. Alignment of vorticity with the eigenvectors of the rate of strain: conditional statistics 

The instantaneous alignment of vorticity with the intermediate eigenvalue of the rate of strain, documented 
many times before for the true velocity gradient tensor m, and for the perceived velocity gradient M (see 
Fig- E}, results from an averaging over all configurations in the flow, which gives the same statistical weight 
to any configuration. This may lead to the impression that, in the case of the velocity gradient tensor, the 
direction of the intermediate eigenvalue of the rate of strain is more important than the largest eigenvalue. 
However, such casual arguments are deceiving^ 6 -, since the comparatively rare configurations where vorticity 
is aligned with the largest eigenvalue of the rate of strain may actually contribute more to the dynamics 
than the frequent contributions, where vorticity is aligned with the intermediate eigenvector. Here we check 
this argument more quantitatively by studying the alignment weighted by the magnitudes of the rate of 
strain and vorticity. Again, we study the evolution in the coordinates of ej(0), the eigenframe of the rate of 
strain at some earlier time. 

Figure [S] presents ([e;(0) ■ u)(t)] 2 ) / (tu 2 {t)), i.e., the evolution of the vorticity weighted alignment with ej(0) 
(or the relative components of enstrophy in the ej(0) direction). Compared with Fig. [6^,, the contribution 
in ei(0) direction is even larger, a consequence of the increase in vorticity due to vortex-stretching in that 
direction^. The contribution in e2(0), on the other hand, decreases with time, which is to be compared 
with the nearly constant ([e2(0) • e u {t)] 2 ) for t/to < 0.2 in Fig. [Bb. 

The alignment weighted by the magnitude of eigenvalues of the initial rate of strain is shown in Fig. [9] 
which calls for the following remarks: (i) The instantaneous alignment is uncorrclated with the magnitude of 
the rate of strain as (Ai(0)[e;(0)-e u (0)] 2 ) w (A i (0))([e. i (0) -e^O)] 2 ). (ii) The increase of (A, (0)[e;(0) -e w (£)] 2 ) 
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FIG. 8. (Color online) Alignment of e„(f) with ei(0) weighted by vorticity. The averages ([ei(0) ■ Lu(t)] 2 )/{uj 2 (t)) are 
shown as a function of t/to for different values of ro in the inertial range. The DNS data correspond to R\ = 170, 
and the experiments to R\ — 350. 



with time in the ei(0) direction appears to be uncorrelatcd with the magnitude of Ai(0) as the same 
amount of increase (approximately 30%) is achieved in the unconditional statistics ([ei(0) ■ e^t)] 2 ) (see 
Fig. [5Ji). Similarly, the contribution in e 2 (0) direction is nearly constant at t/to < 0.2 before decreases, as 
([e 2 (0) • e u {t)} 2 ) in Fig.[6p. 

The quantity (Ai(0)[e;(0)-u;(i)] 2 ) measures the evolution of the enstrophy component in the ei(0) direction, 
weighted by the eigenvalues of the rate of strain, Xi(0). At t = (instantaneous statistics), the quantity 
(Ai(0)[ei(0) • w(0)] 2 ) reduces to the contribution of the i th eigendirection of S to the vorticity growth: 
^i =1 (A i (0)[e. i (0) • w(0)] 2 ) = (uj(0) ■ S • w(0)), which is the usual vortex stretching. Fig. [TUlthus shows the 
contribution to vortex-stretching from the directions of ei(0) and e2(0). It is clear that the contribution 
from ei(0) direction increases rapidly with time and is much larger than that from e2(0). Moreover, the 
contribution in ei(0), once increased, stays at the larger value for longer time as the statistics conditioned 
on vorticity or the rate of strain alone (see Figs.EK and UK), which indicates that larger rate of strain Ai(0) 
results in stronger vorticity increase and persists longer. 
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FIG. 9. (Color online) Alignment of e w (i) withej(O) weighted by strain. The averages (A i (0)[e i (0)-e UJ (t)] 2 >/(aj 2 (0)) 1/2 
are shown as a function of t/to for different values of ro in the inertial range. The DNS data correspond to R\ = 170, 
and the experiments to R\ — 350. 



D. Evolution of the intermediate eigenvalue of the rate of strain 

In addition to the alignment between vorticity and the eigen-directions of the rate of strain, e^, M is also 
characterized by the eigenvalues of the rate of strain, in particular by the intermediate one, A2. Here, we 
consider the normalized value /?, defined by Eq. [13] Fig. [TT] shows the mean value of (3 as a function of 
time t/to, f° r several values of the tetrahedron size, r . Remarkably, the curves superpose extremely well 
for t/to 0.25, as if the evolution were self-similar. 



VII. DISCUSSION AND CONCLUSIONS 



It is appropriate to stress here that asking questions about M is merely one way of studying the correlation 
function of the velocity field u(x, £), based on four different spatial points. Although the precise relation 
between M and the true velocity gradient tensor, m, or its coarse-grained generalization is not completely 
obviousAS, it is our contention that the correlation functions investigated here, based on M, contain far 
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FIG. 10. (Color online) Alignment of e w (t) with ei(0) weighted by both strain and vorticity. The averages (Aj (0) [e* (0) • 
uj(t)] 2 ) / {uj 2 (t)) (^(O) 2 } 1 ^ 2 are shown as a function of t/to for different values of ro in the inertial range. The DNS 
data correspond to R\ — 170, and the experiments to R\ — 350. 



more information than the correlation function of the velocity at two spatial points, such as the structure 
function, and much can be learned about the turbulent flow by investigating correlation functions at more 
than two spatial point o 47 ' 48 . 



A. Comparison between dissipative and inertial scale dynamics 

In the previous sections, the properties of M have been studied for tetrahedra of size ro in the inertial 
range. It is appropriate to compare these alignment properties with those obtained from the true velocity 
gradient tensor, m. We first note that M reduces to the true velocity gradient tensor when r becomes 
significantly smaller than the Kolmogorov scale r\. Here we report results from the DNS on the evolution 
of the alignment between the direction of the true vorticity e w (t) with the eigenvectors ei(0) of s. The 
statistics characterizing the alignment are shown in Fig. [T2l The curves in Fig. [12] and Fig. |6] show strong 
qualitative similarities. They differ only in the scaling of the i-axis. We find that in the dissipative range 
the characteristic time scale for alignment between vorticity and strain is ^ 2 — 3tk , which coincides with 
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FIG. 11. (Color online) Evolution of (/3), the mean of the normalized intermediate eigenvalue of strain. With time 
rescaled with to, the curves corresponding to different values of ro collapse very well. The DNS data correspond to 
Rx — 170, and the experiments to R\ = 350. 



the correlation time of the rate of strain tenso r 49 i 50 , where tj< = ivje) 1 ! 2 is the Kolmogorov time scale. 
Remarkably, the strongest alignment of e w (t) with ei(0) in the dissipative range is as strong as it is in 
the inertial one: compare Fig. and Fig. The tendency of e u (0) to align with e 2 (0), and to be 

perpendicular to 63(0) is stronger in the dissipative range than in the inertial range, as shown also in 
Fig-EJc). Thus, the properties of alignment between e tlJ (t) and ei(0) for M when ro is in the inertial range 
are quite similar with the properties of the true velocity gradient tensor m. Similar observations were 
reported in another very recent numerical study 4 ^. 



B. Fixed shape tetrahedra 

The coupling between the dynamics of M and the geometry of the set of points has been postulated to be 
a crucial ingredient in any effort to properly understand the dynamics of the velocity gradient tenso r 13 ' 18 ' 42 . 
The recent experimental and numerical work 1 ^ effectively demonstrates that the properties of alignment 
of vorticity with the eigenvalues of the rate of strain are strongly dependent on the deformation of the 
tetrahedra. In this context it is worth to ask how the dynamics of alignment is modified, if instead of 
following Lagrangian tetrahedra that are deformed by the flow, one considers tetrahedra of fixed shape and 
size. To analyze this we followed numerically one single tracer particle, and used it as the center of mass of 
an isotropic tetrahedron with a fixed size and arbitrary orientation. A similar approach was used before^ 1 - 
for fixed volumes following Lagrangian particles. For values of ro much smaller than the Kolmogorov length 
scale, 77, the velocity gradient tensor obtained using this construction reduces to the true velocity gradient 
tensor m, whose properties are illustrated in Fig. H121 

Fig. IT51 shows the evolution of the alignment of e^(i) with the eigenvectors of the rate of strain ej(0) for 
the velocity gradient tensor obtained from fixed shape tetrahedra. Comparing Fig. [B] and Fig. Q2] reveals 
some strong quantitative differences. The tendency e u (i) to align with ei(0), as measured by the increase 
of ((e w (i) • ei(0)) 2 ), is much weaker for fixed-shape tetrahedra than for deformable tetrahedra. Similarly, 
the tendency of e w (i) to become more perpendicular to 63(0) is much reduced. We also notice that the time 
scale of the alignment process is much longer for fixed-shape tetrahedra. Moreover, the degrees of alignment 
for fixed-shape tetrahedra (Fig. [T5)) are much weaker than for the dissipative range. 
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FIG. 12. (Color online) Alignment between e w (t) and ei(0) for the true velocity gradient tensor m, obtained from 
DNS (i?A = 170). The averages ([e»(0) ■ e u (t)] 2 } for i = 1, 2 and 3, are indicated by crosses, downwards pointing, 
and right pointing triangles. The time is normalized by the Kolmogorov time scale, tk = (u/e) 1 ^ 2 . The observed 
evolutions of the alignment of e„(£) and ei(0) are qualitatively very similar to what is shown in Fig. [5] for tetrahedra 
with size ro in the inertial range. 



The results shown in Fig. [T3] thus demonstrate a very important feature: the properties of alignment of 
M, constructed with properly deformablc tetrahedra are much closer to those of the true velocity gradient 
tensor, than those obtained with fixed shape tetrahedra. The coupling between the dynamics of M and of 
the geometry is crucial to properly understand and model fluid motion. 

C. Summary 

To sum up, we have investigated systematically the perceived velocity gradient M from the point of view 
of alignment between vorticity and the rate of strain, using both DNS at R\ = 170, over the entire range of 
scales, and experiments at a higher Reynolds numbers (R\ = 350) but over a more restricted range of scales 
in the inertial range. Our results generalize our recent worki^, which focused on the alignment between 
vorticity and the eigenvector ei corresponding to the largest stretching. 

We found that strong deformation of tetrahedra, leading to the formation of flattened, almost coplanar 
configurations, occurs over a time scale of the order to/ 4, which is comparable to the time scale characterizing 
the alignment of vorticity with the eigenvector ei(0) of the strain, corresponding to the largest eigenvalue, 
~ to/5, see Fig. Uta). We note in this respect that during this time, the relative orientation of the vectors 
e u (t) and e 2 (0) remains nearly constant, see Figs.^b) and[7fb), (e), and (h). 

Whereas the dynamics leading to alignment between uj and ei has been found to be essentially self- 
similar— over the whole range of scales and Reynolds numbers considered, the self-similarity of the statistical 
properties of the angle between uj and the other two eigendircctions of the rate of strain is only observed at 
the largest Reynolds number studied here, R\ = 350. At small Reynolds numbers, R\ < 170, we found a 
systematic dependence on the size ro of the tetrahedra of the alignment of vorticity and the rate of strain 
at the same time. We also observed that the evolution of the angle between vorticity and the second and 
third eigenvalues of the rate of strain does not simply reduce to a simple function of t/to- In any event, the 
dynamics proceed with a characteristic time scale of order ~ to/ 5. 

For the range of Reynolds numbers considered here, the picture that emerges from the study of M 
reflects the properties in the dissipative scale in terms of alignment of vorticity and strain eigenvectors. The 
characteristic time of evolution in the dissipative scale is ~ 2 — 3tk ■ 

Future studies will be devoted to a more systematic comparison with modeling, both in terms of the 
tetrahedron models, and in terms of the arguments developed in our recent worki^. 
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FIG. 13. (Color online) Alignment between e u (t) and ei(0) measured with isotropic tetrahedra of fixed shapes and 
size ro, obtained from DNS at R\ = 170. The averages <[ei(0) ■ e„(t)] 2 ) for i = 1 (a), 2 (b) and 3 (c), as a function 
of t/to, are shown for several values of ro, all in the inertial range. The tendency of e^(t) to align with ei(0), and to 
become perpendicular to e3(0), are much weaker than with tetrahedra freely advected with the flow (compare with 
Fig. O, or with the true velocity gradient tensor (compare with Fig. 1121) . 
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